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1. Introduction 

Supersymmetry (SUSY) has been the top eandidate for beyond standard model (BSM) physies 
sinee many years. This has several reasons: SUSY solves the hierarchy problem of the standard 
model (SM), provides a dark matter candidate, leads to gauge coupling unification, and gives an 
explanation for electroweak symmetry breaking (EWSB). Before the LHC has turned on, the main 
focus has been on the minimal supersymmetric extensions of the SM, the MSSM. However, the 
negative results from SUSY searches at the LHC as well as the measured Higgs mass of about 
125 GeV put large pressure on the simplest scenarios. Wide regions of the parameter space, which 
had been considered as natural before LHC has started, have been ruled out. This has caused 
more interest in non-minimal SUSY models. Also non-SUSY extensions of the SM became more 
popular again. Other BSM models can provide many advantages compared the MSSM: 

• Naturalness: about one third of the Higgs mass has to be generated radiatively in the MSSM 
to explain the observation. Heavy SUSY masses are needed to gain sufficiently large loop 
corrections, i.e. a soft version of the hierarchy problem appears again. The need for large 
loop corrections gets significantly softened if F- or D-terms are present which already give 
a push to the tree-level mass. 

• SUSY searches: the negative results from all SUSY searches at the LHC have put impressive 
limits on the Sparticle masses. However, the different searches are based on certain assump¬ 
tions. As soon as these conditions are no longer given like in models with broken /^-parity, 
the limits become much weaker. 

• Neutrino data: neutrino masses are not incorporated in the MSSM. To do that, either one 
of the different seesaw mechanisms can be utilised or /^-parity must be broken to allow a 
neutrino-neutralino mixing. 

• Strong CP-problem: the strong CP problems remains not only an open question in the SM 
but also in the MSSM. In principle, for both models the same solution exists to explain the 
smallness of the 0 term in QCD: the presence of a broken Peccei-Quinn symmetry. 

• /t -problem: the superpotential of the MSSM involves one parameter with dimension mass: 
the /r-term. This term is not protected by any symmetry, i.e. the natural values would be 
either exactly 0 or 0{Mgut)- However, for phenomenological reasons it should be similar 
to the electroweak scale. This could be obtained if the /r-term is actually not a fundamental 
parameter but generated dynamically like in singlet extensions. . 

• UV completion: it is not clear that only the gauge sector and particle content of the MSSM 
is the low energy limit of a GUT or String theory. Realistic UV completions predict often 
many additional matter at the TeV scale. In many cases also additional neutral and even 
charged gauge bosons are present. 

• R-symmetry: if one considers R-symmetric models, Majorana masses for gauginos are for¬ 
bidden. To give masses to the gauginos in these models, a coupling to a chiral superfield 
in fhe adjoinf represenfafion is needed. This gives rise fo Dirac masses for fhe gauginos 
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which are in agreement with /^-symmetry. Dirac gauginos are also attractive because they 
can weaken LHC search bounds, and flavour constraints. 

Despite the large variety and flexibility of SUSY, many dedicated public computer tools like 
SoftSUSY, SPheno, Suspect, Isa jet or FeynHiggs are restricted to the simplest realiza¬ 
tion of SUSY, the MSSM, or small extensions of it. Therefore, more generic tools are needed to 
provide the possibility to study many other models with the same precision as the MSSM. This pre¬ 
cision is needed to confront also these models with the strong limits from collider searches, flavour 
observables, dark matter observations, and Higgs measurements. Such a tool is the Mathematica 
package SARAH. I’ll given in the next section an introduction to SARAH, before I explain the in¬ 
terface to SPheno in sec. 3, to Monte-Carlo (MC) tools in sec. 4, and to other tools in sec. 5. I 
summarize in sec. 6. Interested readers find many more details about the presented framework in 
Ref. [1]. 

2. BSM models and SARAH 

2.1 Model definition 

2.1.1 Input needed by SARAH to define a model 

SARAH [2, 3, 4, 5, 6] is optimized for the handling of a wide range of BSM models. The 
basic aim of SARAH is to give the user the possibility to implement models in an easy, compact 
and straightforward way. Most tasks to get the Lagrangian are fully automatized; it is sufficient to 
define just the fundamental properties of the model: 

1. Global symmetries 

2. Gauge symmetries 

3. Chiral superfields 

4. (Super)potential 

5. Field rotations 

That means that SARAH automatizes many steps to derive the Lagrangian from that input: 

1. All interactions of matter fermions and the F-terms are derived from the (superjpotential 

2. All vector boson and gaugino interactions as well as D-terms are derived from gauge invari¬ 
ance 

3. All gauge fixing terms are derived by demanding that scalar-vector mixing vanishes in the 
kinetic terms 

4. All ghost interactions are derived from the gauge fixing terms 

5. All soft-breaking masses for scalars and gauginos as well as the soft-breaking counterparts 
to the superpotential couplings are added automatically for SUSY models 
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6. All rotations from the gauge to the mass eigenstates are performed 
The range of models which can be handled by SARAH is very broad: 

• Global symmetries: SARAH can handle an arbitrary number of global symmetries which 
are either or [/(I) symmetries. Also a continuous 7?-symmetry f/(l)s is possible. 

• Gauge sector: SARAH is not restricted to the SM gauge sector, but many more gauge groups 
can be defined. To improve the power in dealing with gauge groups, SARAH has linked 
routines from the Mathematica package Susyno [7]. SARAH together with Susyno 
take care of all group-theoretical calculations: the Dynkin and Casimir invariants are cal¬ 
culated, and the needed representation matrices as well as Clebsch-Gordan coefficients are 
derived. Also gauge kinetic mixing is fully supported for an arbitrary number of Abelian 
gauge groups. 

• Matter sector: there is no practical restriction how many matter states -either defined as chi¬ 
ral superfields in SUSY models, are as component fields in non-SUSY models- are present. 
All states can come with an arbitrary number of generations and can transform as any irre¬ 
ducible representation with respect to the defined gauge groups. 

• (Super)potential: the main restriction is that only renormalizable terms are supported in the 
(super)potential. However, there is no restriction how many terms are present. 

As example how short a model implementation in SARAH is, I show in Fig. 1 the model file for the 
MSSM. 

2.1.2 Checks of implemented models 

SARAH provides functions to check the (self-) consistency of a model addressing the following 
questions: 

• Causes the particle content gauge anomalies? 

• Leads the particle content to the Witten anomalies? 

• Are all terms in the (super)potential in agreement with global and local symmetries? 

• Are there other terms allowed in the (super)potential by global and local symmetries? 

• Are all unbroken gauge groups respected? 

• Are there terms in the Lagrangian of the mass eigenstates which can cause additional mixing 
between fields? 

• Are all mass matrices irreducible? 

• Are the properties of all particles and parameters defined correctly? 
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Off{GeneraI: :$peUj 

MtodelNam = 'HSSM'; 

Moder NaneiaTeX ='IISSM’; 
Model Authors - *F.Staub’; 
HocterDate * 


(• Particle Content*) 

(•.•) 

(* Gauge Superfields *) 

GaugeC(']]-{B. U[t], hypercharge, g1.False); 
6auge[[;i3»{«. SUt2]. left. g2.True): 

Gauged.'DMC. SU[)], color. g3,False): 


(* Oiiral Superfields *) 

FieldstO) » ffuL. dL), J. q, 1/6. 2. 3) 
Fieldstt.'ll » {M, eL), 1- -U2. 2. 1> 
Fields[[?]] = {{HdO. Hdn). , Hd. -1/2. 2, 1) 
Fields{[J]] = {{Hup. HuO), Hu. 1/2. 2. ') 


Fields{[EJ) » {conjEdR]. 1. d. VS, -3): 

Fields[(i3) * {conjtuRJ. 3. u. -2/3. -3): 

Fjelds{(^l) « {conj[eRI. 3 . e. '); 

(•.•) 

(• Superpotential •) 

(•.•) 

SuperPotential ={{{'. Yu),{u,q.Hu)). {{-•.Yd).{d.q.Hd)), 
{{-•.Ye).{e,l.Hd)). {{■.\(llu]).{Hu.Hd))] 

(•.•) 

(* Integrate Out or Delete Particles *) 

(•.•) 


lntegrateOot*{); 

OeleteParticles-{): 


NaMOfStates={GaugeES. EWSfi), 


(•.After EWSB.•) 


(• Gauge Sector •) 

DEFINITIONtEWSeKGauieSectorl « 

{ 

{{VB.V«[-'I),{VP.V2}.ZZ). 

{{Wfi [ • 1. VHB { : 3). {Vkpi. coo j [ W*I ]), ar). 

{{ fHB [ ‘}. fMB {.]. ne [:]}.{ fMe, fWp, fM)). Z fM) 

); 


(».VEVs •} 

D€FINITIOIC[EWSB] {VEVs]= 

{{SHdO. {vd, '/Sqrt[4). {sigead, \[lMginaryI]/Sqrti:]).{phid.i/Sqrt[2])). 
{SHoO. {vu. '/Sqrttii). {sigeau, vtI»agifkaryli/Sqrtt:j).{phiu.'/Sqrtt2j))): 


(•-Kixings-•) 

0€FINITI{M[EWSB](llatterSector]: 

{ {{SdL. SdR). <Sd, ZO)). 

|{SvL). {Sv. ZV)). 

{{SuL. SuR). {Su. ZU)). 

{{SeL. SeR). {Se. ZE)). 

{{phid. phiu). {hh. ZH)), 

{{signad. signau). {Ah. ZA)). 

{{SHd»,conj{SHupl).{Hp«.ZP)). 

{{fB. «0, FHdO, FHuO), {LO. ZN)>. 

{{{fV*i. FHd*). {flep, FHup)). {{Lji.UII). {Lp.UP))). 
{{{FeL).{conj[FeRI)),{{FEL,ZEL).{FER.ZER))). 
{{{FdL).{C0fij[FdR])).{{Ft)t,2DL).{FDR.Z»))). 
{{{Fut).{coojIFuR])).<{FUl.ZUL).{F«l.ZUR))) 

>; 

DEFINITIQNIEMSB]tPhases]= 

{ {fG. PhaseGlu) ); 


OCFINITIONIEMSB] {OiracSpir>ors]s{ 
Fd ->{ FDl. coojIFORJ), 

Fe ->{ FEL, conjlFERJ), 

Fu ->{ FUL. conjlFURJ), 

Fw ->{ FvL. :), 

Chi ->< LO, conjtLOJ). 

Cha •>< La. conj[Lpi). 

Glu ->{ fG. conj[fG]) 

); 


Figure 1: The model file to implement the MSSM in SARAH 

2.2 Calculations and output 

SARAH can perform in its natural Mathematica environment many calculations for a model 
on the analytical level. For an exhaustive numerical analysis usually one of the dedicated interfaces 
to other tools discussed in the next sections is the best approach. I start with some details about the 
analytical calculations which SARAH can perform. 

2.2.1 Tree-level properties 

Tadpole equations During the evaluation of a model, SARAH calculates ’on the fly’ all minimum 
conditions of the tree-level potential, the so called tadpole equations. 

Masses SARAH calculates for all states which are rotated to mass eigenstates the mass matrices 
during the evaluation of a model. In addition, it calculates also the masses for the states for which 
no field rotation has taken place. 

Vertices SARAH has functions to extract in an efficient way all tree-level vertices from the La- 
grangian. These vertices are saved in different Mathematica arrays according to their generic 
type. 

2.2.2 Renormalization group equations 

SARAH calculates for SUSY and non-SUSY models the full two-loop RGBs including the full 
CP and flavour structure. I give here a short summary of the underlying, generic results which are 
used in these calculations. 
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SUSY RGEs The calculation of the SUSY RGBs is mainly based on Ref. [8]. However, this 
work does not cover all possible subtleties which can appear in SUSY models. Therefore, SARAH 
has implemented also some more results from literature which became available in the last few 
years: 


• In the case of several 17(l)’s, gauge-kinetic mixing can arise if the groups are not orthogonal. 
Substitution rules to translate the results of Ref. [8] to those including gauge kinetic mixing 
where presented in Ref. [9] and are used by SARAH 

• The calculation of the RGEs in the presence of Dirac gaugino is based on Ref. [10] 

• The results of Refs. [11, 12] are used to get the gauge dependence in the running of the 
VEVs. 

Non-SUSY RGEs SARAH uses the expressions of Refs. [13, 14, 15, 16] for the calculation of the 
RGEs in a general quantum field theory. These results are completed by Ref. [17] to cover gauge 
kinetic mixing and again by Refs. [11, 12] to include the gauge-dependence of the running VEVs 
also in the non-SUSY case. 

2.2.3 One- and two-loop corrections to tadpoles and self-energies 

One-loop corrections SARAH calculates the analytical expressions for the one-loop corrections 
to the tadpoles and the one-loop self-energies for all particles. Eor states which are a mixture of 
several gauge eigenstates, the self-energy matrices are calculated. The calculations are performed 
in DR-scheme using’t Hooft gauge. In the case of non-SUSY models SARAH switches to MS- 
scheme. All generic diagrams which are covered at one-loop are depicted in Eig. 2. 



Figure 2: Generic diagrams included by SARAH to calculate one-loop tadpoles and self-energies. 


Two-loop corrections It is even possible to go beyond one-loop with SARAH and to calculate 
two-loop contributions to the self-energies of real scalars. There are two equivalent approaches 
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implemented in the SPheno interface of SARAH to perform these calculations: an effective poten¬ 
tial approach, and a diagrammatic approach with vanishing external momenta. More details about 
these calculations are given in sec. 3.1.3. 

2.3 Output 

The information derived by SARAH can be outputted in different ways 

• Model files for Monte-Carlo tools like CalcHep, MadGraph, or WHIZARD 

• Model files for ofher fools like FeynArts and Vevacious 

• Forfran code fo implemenf a model in SPheno^. 

• FTpX files 

A few more defails abouf fhe differenl oufpufs are given in fhe nexf three sections. 

3. Linking SPheno and SARAH 

SARAH 3 has been the first ‘spectrum-generator-generator’: SARAH writes Fortran source 
code for SPheno using the derived information about the mass matrices, tadpole equations, ver¬ 
tices, loop corrections and RGBs for the given model. Based on this code the user gets a fully 
functional spectrum generator for a model. The features of a spectrum generator created in this 
way are 

• Full two-loop running of all parameters and all masses at one-loop 

• two-loop corrections to Higgs masses 

• Complete one-loop thresholds at Mz 

• calculation of flavour and precision observables af full 1-loop 

• calculafion of decay widfhs and branching rafios 

• interface fo HiggsBounds and HiggsSignals 

• esfimafe of elecfroweak Fine-Tuning 

I give a few more defails abouf fhe differenf calculations. 

3.1 Mass calculation with SPheno 

3.1.1 Threshold corrections 

For a precise calculation of the masses, it is necessary to have an accurate input for all param¬ 
eters which enter the calculation. In general, the running SM parameters depend on the masses of 
BSM states. The reason are the thresholds to match the running parameters to the measured ones. 
The routines generated by SARAH perform a full one-loop matching in the given model to calculate 
the SM gauge and Yukawa couplings. 

• A similar interface to generate a spectrum generator written in C+ became available via FlexibleSUSY [18]. 


7 




Introduction to SARAH and related tools 


Florian Staub 


3.1.2 One-loop masses 

The one-loop mass spectrum is calculated from the information about the one-loop self¬ 
energies and tadpole equations. The procedure is a generalisation of the one explained in detail 
for the MSSM in Ref. [21]. The main features are 

• Any one-loop contribution in a given model to all fermions, scalars and vector bosons is 
included 

• The full dependence in the loop integrals is included 

• An iterative procedure is applied to find the on-shell masses m{p^ = m^). 

3.1.3 Two-loop Higgs masses 

SARAH can also generate Fortran code to calculate the two-loop corrections to real scalars 
with SPheno. The same approximation as usual done for the MSSM are also applied here: (i) all 
calculations are performed in the gaugeless limit, i.e. the electroweak contributions are dropped, 
and (ii) the momentum dependence is neglected. Using these routines,the theoretical uncertainty in 
the Higgs mass prediction for many models has shrunken now to the level of the MSSM. In general, 
there are two different techniques to calculate the two-loop corrections with SARAHSPheno: 

• Effective potential calculation: SARAH makes use of the generic two-loop results for the 
effective potential given in Ref. [22]. To get the values for the two-loop self-energies and 
two-loop tadpoles, the derivatives of the potential with respect to the VEVs are taken numer¬ 
ically as proposed in Ref. [23]. There are two possibilities for this derivation implemented in 
SARAH/SPheno: (i) a fully numerical procedure which takes the derivative of the full effec¬ 
tive potential with respect to the VEVs. A semi-analytical derivation which takes analytical 
the derivative of the loop functions with respect to involved masses, but derives the masses 
and coupling numerically with respect to the VEVs. More details about both methods and 
the numerical differences are given in Ref. [24]. 

• Diagrammatic calculation: A fully diagrammatic calculation for two-loop contributions to 
scalar self-energies with SARAH-SPheno became available with Ref. [25]. In this setup a 
set of generic expressions first derived in Ref. [25] is used. All two-loop diagrams shown in 
Eig. 3 are included in the limit p^ = 0. The advantage of the diagrammatic approach is that 
no numerical derivation is needed. 

A comparison between the results using the new routines generated by SARAH and widely used 
routines based on Refs. [26, 27, 28, 29, 30] for the MSSM and Ref. [31] for the NMSSM is shown 
in Eig. 4. Obviously, a perfect agreement is found for both models. Since there are no further 
possibilities to cross check the results obtained for the Higgs mass with SARAH/SPheno, the above 
described implementation of two independent approaches (effective potential, diagrammatic) was 
necessary to have the possibility to double check results. In the meantime, this setup has been used 
to analyse the Higgs mass also in other non-minimal SUSY models at the two-loop level: 

• Contributions of trilinear RpV, see Ref. [32] 



Introduction to SARAH and related tools 


Florian Staub 



Figure 3: Generic diagrams included to calculate the two two-loop corrections to tadpoles and scalars. 
These are the diagram which don’t vanish in the gaugeless. 


• Missing corrections for NMSSM, see Ref. [33] 

• Contributions from non-holomorphic soft-terms, see Ref. [34] 

• MRSSM at full two-loop, see Ref. [35] 

• Contributions from vectorlike stops, see Ref. [36] 

3.2 Decay widths and branching ratios 

SPheno modules created by SARAH calculate all two-body decays for SUSY and Higgs states 
as well as for additional gauge bosons. In addition, the three-body decays of a fermion into three 
other fermions, and of a scalar into another scalar and two fermions are included. 

In the Higgs sector, possible decays into two SUSY particles, leptons and massive gauge bosons 
are calculated at tree-level. For two quarks in the final state the dominant QCD corrections due to 
gluons are included [37]. The loop induced decays into two photons and gluons are fully calculated 
at leading-order (LO) with the dominant next-to-leading-order corrections known from the MSSM. 
For the LO contributions all charged and coloured states in the given model are included in the 
loop, i.e. new contributions rising in a model beyond the MSSM are fully covered at one-loop. In 
addition, in the Higgs decays also final states with off-shell gauge bosons (ZZ*, WW*) are included. 
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Figure 4: Comparison of the Higgs mass calculation at two-loop level using the automatically generated 
routines of SARAH (full lines) and the routines based in Refs. [26, 27, 28, 29, 30] for the MSSM and Ref. [31] 
for the NMSSM. The green lines include only the corrections 0{as{a, + «/>)), while for the blue lines all 
other, available corrections are included. The black line corresponds to the one-loop results. 

3.3 Flavour observables 

The SPheno modules written by SARAH contain already out of the box the routines to calcu¬ 
late many quark and lepton flavour violating observables: 

• Lepton flavour violation: 

- Br(/i ^ IjY), Br(/ ^ 3/'), Br(Z ^ W) 

- CR{n-e,N) (N=Al,Ti,Sr,Sb,Au,Pb),Br(T ^Z + P) (P=7i, T],T]') 

• Quark flavour violation: 

- Br(P ^ X,Y), Br(p0^ ^ ZZ”), Br(P ^ slf), Br(K juv) 

- Br(B —^vv), Br(isr+ —^ tt+vv), Br(Ki — tt^vv) 

- AMb,,b,, AM/f, Ek, Br(P ^ KnfL) 

- Br(B ^ Zv), Br(D, ^ Zv) 

The calculation is based on the FlavorKit functionality [38] which is a connection FeynArt s- 
SARAH-SPheno. This provides a full one-loop calculation in a given model. In addition, this 
interface can be used to derive Wilson coefficients for new operators, and to relate coefficients 
to calculate new observables with SPheno. In order to demonstrate that this fully automatized 
calculation reproduces the results from dedicated MSSM codes, a comparison is shown in Fig. 5. 


10 











Introduction to SARAH and related tools 


Florian Staub 




tan/3 


mo [GeV] 




Afi /2 [GeV] 


3\/i/ 2 [GeV] 


Figure 5: Comparison of different flavour observables in the CMSSM. The colour code is: FlavorKit 
(red), SPheno 3.3 (blue), SUSY_Flavor 1 (brown), SUSY_Flavor 2 (green), MicrOmegas (or¬ 
ange), Super I so (plum). 

3.4 Checking Higgs constraints 

HiggsBounds [39,40,41] and HiggsSignals [42] are dedicated tools to study the Higgs 
properties of a given parameter point in a particular model. While HiggsBounds checks the pa¬ 
rameter point against exclusion limits from Higgs searches, HiggsSignals gives a value to 
express how good the point reproduces the Higgs measurements. In general, HiggsSignals 
and HiggsBounds can handle different inputs: either the cross sections for all necessary pro¬ 
cesses can be given at the parton or hadron level, or the effective couplings of the Higgs states 
to all SM particles are taken as input. In addition, the masses and widths of all neutral as well 
as charged Higgs states are always needed. SPheno provides all input for the effective coupling 
approach. The information is given in the SLHA spectrum file and in addition in separated files 
(called MH_GammaTot. dat, MHplus_GammaTot. dat, BR_H_NP . dat, BR_Hplus . dat, 
BR_t. dat, ef fC . dat). While SLHA files can be used with HiggsBounds for models with 
up to five neutral scalars, the separated files can be used with even up to nine neutral and nine 
charged scalars. 

4. MC Tools and SARAH 

SARAH writes all necessary files to implement a new model in different MC tools. I give a 
short overview of the main features and explain, how parameter values between SPheno and the 
MC can be exchanged. 
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4.1 Output of models files 

CalcHep The model files for CalcHep [43, 44] can also be used with MicrOmegas [45] 
for dark matter studies. CalcHep is able to read the numerical values of the masses and mixing 
matrices from a SLHA spectrum file based on the SLHA+ functionality [46]. SARAH makes use of 
this functionality to generate the model files in a way that they read the spectrum files written by a 
SARAH generated SPheno version. 

UFO format SARAH can write model files in the UFO format [47] to implement new models 
in MadGraph 5. This format is also supported by other tools like GoSam [48], Aloha [49], 
Herwig++ [50, 51] or MadAnalysis 5 [52]. 

One can give the spectrum file written by SPheno as parameter card to MadGraph. 

WHIZARD/O'Mega SARAH writes all necessary files to implement a model in WHIZARD and 
O' Mega [53, 54]. Since the SLHA reader of WHIZARD is at the moment restricted to the MSSM 
and NMSSM, SPheno versions generated by SARAH write all information about the spectrum and 
parameters also in an additional tile in the WHIZARD specific format. This tile can be then read by 

WHIZARD. 

4.2 Interplay SARAH-SPheno-MC-Tool 

The tool chain SARAH-SPheno-MC-Tools has one very appealing feature: the implementa¬ 
tion of a model in the spectrum generator (SPheno) as well as in a MC tool is based on just one 
single implementation of the model in SARAH. Thus, the user has not to worry that the codes might 
use different conventions to define the model. In addition, SPheno also provides all widths for the 
particles and this information can be used by the MC-Tool as well to save time. 

5. Other output 

5.1 FeynArts/FormCalc 

FeynArts/FormCalc [55, 56] is a very powerful combination of two Mathematica 
packages to perform tree-level and loop calculations. FeynArts is a Mathematica package for 
the generation and visualization of Feynman diagrams and amplitudes. FormCalc makes use of 
Feynman diagrams generated by FeynArts, and of Form for analytical manipulations: it calcu¬ 
lates the amplitude, squares it, performs colour simplifications, makes the polarisation sums, and 
so on. FormCalc can also write Fortran or C code for further numerical calculations. SARAH 
can be used to generate the model files for these tools to implement a new model. 

5.2 Vevacious 

Vevaclous is a tool to check for the global minimum of the one-loop effective potential for 
a given model allowing for a particular set of non-zero VEVs. For this purpose Vevacious finds 
first all tree-level minima by using HOM4PS2 [57]. Afterwards, it minimizes the one-loop effective 
potential starting from these minima using minuit [58]. If the input minimum turns out not to 
be the global one, life-time of meta-stable vacua can be calculated using Cosmotransitions 
[59]. 
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Vevacious takes the tadpole equations, the polynomial part of the scalar potential and all mass 
matrices as input. All of this information has to be expressed including all VEVs which should 
be tested. That means, in order to check for charge and colour breaking minima the stop and stau 
have to show up in the potential and mass matrices and the entire mixing triggered by these VEVs 
should be included. To take care of all that, the corresponding input files can be generated by 
SARAH. Possible applications when these checks become important are: 

• Charge and colour breaking minima, see for instance Refs. [60, 61] 

• Spontaneous R-parity violation, see for instance Ref. [62] 

• NMSSM with large A terms, see for instance Refs. [63, 64] 

5.3 BTeX 

All derived information (mass matrices, vertices, RGEs, loop corrections) can be exported 
to ETgX to get the expressions in a human readable form. There exists the possibility to draw 
Eeynman diagrams with FeynMF [65] for all vertices. 


6. Summary 


I have given in this lecture a short summary how the study of a new particle physics models 
can be automatized to a large extent. The first step is to implement a new model in SARAH. SARAH 
then calculates all analytical properties of this model and passes this information to other tools for 
the numerical evaluation. The main work-flow of this setup is summarized in Eig. 6. 
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Figure 6: Schematic framework to study all aspects of a new model in a highly automatized way. 
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